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molecules on surfaces 
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We present a modification of the ASCF metfiod of calculating energies of excited states, in order 
to make it applicable to resonance calculations of molecules adsorbed on metal surfaces, where the 
molecular orbitals are highly hybridized. The ASCF approximation is a density functional method 
closely resembling standard density functional theory (DFT), the only difference being that in ASCF 
one or more electrons are placed in higher lying Kohn-Sham orbitals, instead of placing all electrons 
in the lowest possible orbitals as one does when calculating the ground state energy within standard 
DFT. We extend the ASCF method by allowing excited electrons to occupy orbitals which are linear 
combinations of Kohn-Sham orbitals. With this extra freedom it is possible to place charge locally 
on adsorbed molecules in the calculations, such that resonance energies can be estimated, which 
is not possible in traditional ASCF because of very delocalized Kohn-Sham orbitals. The method 
is applied to N2, CO and NO adsorbed on different metallic surfaces and compared to ordinary 
ASCF without our modification, spatially constrained DFT and inverse-photoemission spectroscopy 
(IPES) measurements. This comparison shows that the modified ASCF method gives results in close 
agreement with experiment, significantly closer than the comparable methods. For N2 adsorbed on 
ruthenium (0001) we map out a 2-dimensional part of the potential energy surfaces in the ground 
state and the 27r-resonance. From this we conclude that an electron hitting the resonance can induce 
molecular motion, optimally with 1.5 eV transferred to atomic movement. Finally we present some 
performance test of the ASCF approach on gas- phase N2 and CO, in order to compare the results 
to higher accuracy methods. Here we find that excitation energies are approximated with accuracy 
close to that of time-dependent density functional theory. Especially we see very good agreement 
in the minimum shift of the potential energy surfaces in the excited state compared to the ground 
state. 

PACS numbers: 31.15.xr, 31.50.Df, 82.20.Gk 



I. INTRODUCTION 

Density functional theory^ '^ (DFT) has proved to be a 
vital tool in gaining information on many gas-surface pro- 
cesses. This may be surprising, since DFT is only valid 
for relaxed systems in their ground state and therefore 
not directly applicable to dynamical situations. How- 
ever, often the electrons relax much faster than the time 
scale of the atomic movement, such that the electron gas 
can be considered relaxed in its ground state at all times. 
Then potential energy surfaces (PES) of the ground state 
obtained by DFT, or any other method, can be used to 
describe the motion of atomic cores. This is the Born- 
Oppenheimer approximation. 

In some situations, however, the Born-Oppenheimer 
approximation is not valid. This is for example the 
case when the electronic system is excited by a fem- 
tosecond laser.2.!i. or hot electrons are produced with a 
metal-insulator-metal junctioui^ The Born-Oppenheimer 
approximation also breaks down if the time scales for the 
electronic and nuclear motions are comparable or if the 
separations between the electronic states are very small, 
such that transitions between the electronic states will 
occur. In these situations it is necessary to go beyond 
the Born-Oppenheimer approximation either by consid- 



ering the coupling between electronic states^ where it 
becomes necessary to obtain PESs of excited states, or 
by an electronic friction model i^*^ 

The problem of calculating excitation energies are be- 
ing approached in many different ways, even within DFT. 
Time dependent density functional theory (TDDFT)^- 
gives, compared to the computational cost, good agree- 
ment with experiments for excitations in atoms and 
molecules,— However, TDDFT suffers some problems 
in excitations involving charge transferrin The GW 
approximatio n^^d^ can be used to gain accurate exci- 
tation energies for molecules and clusters. The embed- 
ding method j^^d^ which combines high-accuracy quan- 
tum chemistry methods with DFT, makes it possible to 
handle larger periodic systems with great accuracy. The 
embedding theory has been applied to estimate PESs of 
excited molecules on surfaces.!^ However, the computa- 
tional cost and involved complexity is still very high. Our 
aim has been to find a method, which at a computational 
cost close the level of ground state DFT can estimate 
excited state energies of molecules on surfaces with rea- 
sonable accuracy. Such a method would make it more 
feasible to consider a large range of systems in search of 
systems with interesting or desired properties. 

Constrained DFTlii^ii^ and ASCF^.^ are two differ- 



ent approaches, which both can be considered as small 
extensions of ground state DFT, such that the compu- 
tational cost lies close to that of ground state DFT. In 
constrained DFT an additional potential is introduced 
and varied until a certain constraint on the electrons is 
fulfilled. The simplest approach is to lower (or increase) 
the potential in a certain part of space until you have 
the desired number electrons in this area.^* A different 
approach is to introduce potentials on the orbitals in a 
localized basis set, which depends on the orbitals' po- 
sitions in space.-^ In section Hill we will argue that when 
considering molecular resonance states on surfaces it may 
be problematic with such a strict constraint on the elec- 
trons, since a part of the charge may return to the surface 
on a much shorter time-scale than the lifetime of the res- 
onance. 

In the ASCF scheme the positions of the electrons are 
controlled by controlling the occupation of the Kohn- 
Sham (KS) states as the system reaches self-consistency. 
The ASCF scheme has for a long time been justified in 
cases, where the excited state corresponds to the low- 
est state of a given symmetry. ^2, The scheme has, how- 
ever, often been applied to more general cases. More re- 
cently, G6rling2^ extended the KS formalism to include 
excited states, such that ASCF gets a formal justification 
in the general case, although a special unknown orbital- 
dependent exchange-correlation potential should be used 
for the excited states. In practical implementations stan- 
dard exchange-correlation potentials from ground state 
DFT are typically used. 

This traditional way of just controlling the occupa- 
tion of the KS orbitals has some limitations. For exam- 
ple when a molecule is placed on a metallic surface the 
molecular orbitals will hybridize with the orbitals in the 
surface, such that the molecular orbitals will be spread 
over several KS states. For such systems there is no good 
way of representing a resonance on the molecule as a 
change in the occupations of the KS orbitals. The opti- 
mal thing one can do within this scheme is to occupy the 
KS orbital with the largest overlap with the molecular 
orbital in question, but this overlap can be quite small 
and highly system size dependent. This problem was also 
pointed out by Hellman et al.?''^ and Behler et al^ 

In this paper we modify the ASCF approach, such that 
electrons are allowed to occupy arbitrary linear combi- 
nations of KS orbitals. In this way one achieves much 
better control on the position of the excited electron. As 
is the case for traditional ASCF some knowledge of the 
resonance one wants to consider is needed in order to 
apply the method. The method is especially relevant in 
Newns-Andersson?^i^ type systems, where a resonance 
can be attributed to a known single level, which has been 
hybridized through interactions with other states. This 
includes systems with molecules adsorbed on metal sur- 
faces and molecules trapped between to metal contacts. 

The modification we propose only has minor implica- 
tions on the way practical calculations are performed, 
which is very similar to performing an ordinary ground 



state DFT calculation. In the following we will go 
through the details of the method and apply it to a few 
diatomic molecules on metallic surfaces. The obtained 
results will be compared to the ordinary ASCF method, 
spatially constrained DFT and IPES measurements. Fi- 
nally we present some tests on the performance of the 
ASCF approach on N2 and CO in the gas phase. 



II. METHOD 

In the following we go through the differences between 
the linear expansion ASCF method we propose, ordinary 
ASCF and standard DFT. We start by stating the modi- 
fication of the KS equations when considering an electron 
excited from the Fermi level to a higher lying state. Then 
we show how this affects the energy calculation. Finally 
we generalize the approach to other types of excitations. 



A. Kohn-Sham equations 

The fundamental KS equations^ represent a practical 
way of finding the ground state electron density for a 
given external potential and a given number of electrons 
through an iterative process. 
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where vks is the KS potential, E^c is the exchange- 
correlation energy and N is the number of electrons. As 
seen from Eq. ^ only the N orbitals with lowest energy 
contribute to the density, ie. the electrons are placed in 
these orbitalsi^ In ordinary ASCF one estimates proper- 
ties of excited states by placing the electrons differently. 
For example the HOMO-LUMO gap in a molecule could 
be estimated by replacing Eq. ^ with 
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where V'aCr) is the KS orbital resembling the LUMO from 
the ground-state calculation. Naturally, the KS orbitals 
found when solving these modified KS equations will dif- 
fer from the ones found in an ordinary DFT calculation, 
due to the change in the Hamilton through the change in 
the density when different orbitals are occupied. 

In the linear expansion ASCF method we propose, the 
excited electron is not forced to occupy a KS orbital, but 
can occupy any orbital that is a linear combination of 



empty KS orbitals: 



M 



(5) 



=N 



where M is the number of KS orbitals in the calcula- 
tion. In practice this means that the KS many-particle 
wavefunction is no longer just a Slater determinant of N 
KS orbitals, but a Slater determinant of A^ — 1 KS or- 
bitals and T/;'"'^'* (r) . Only empty KS orbitals are included 
in the linear expansion, since otherwise ^'''^^(r) will not 
be orthogonal to the filled KS orbitals. Eq. ([2]) is then 
replaced with 
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r^{r)Mr)+ J2 a*a,^:ir)i,,{r). (6) 



i,j = N 



Since the expansion coefficients, a^, in principle could 
have any value some a priori knowledge is needed in or- 
der to choose good values. In the case of molecular reso- 
nances on surfaces the expansion coefficients are chosen 
such that ip^'^'^ (r) resembles the relevant molecular orbital 
as much as possible, ie. 



(^^10) 
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where (f> is the molecular orbital. This is consistent with a 
Newns-Andersson2ii^ picture, where the resonance cor- 
responds to an electron getting in the molecular orbital, 
but the resonance broadening and energy shift is due 
to hybridization with the metallic bands and an image 
charge effect. 

In calculations with k-point sampling the linear ex- 
pansion is performed independently in all k-points. In 
the linear expansion ASCF one then avoids the difficul- 
ties one can encounter in choosing which KS state to oc- 
cupy in each k-point in the traditional way of performing 
ASCF calculations: For example one may risk occupying 
different bands in each k-point, when just choosing the 
KS orbital with the largest overlap with the molecular 
orbital in each k-point. 



B. The energy 



nary DFT the kinetic energy is evaluated as 
T[n{v)] = ^(^,1-1^,) 
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where the last equality is seen directly from Eq. ([T]) . Sim- 
ilarly the expression for the kinetic energy in the linear 
expansion ASCF is found to be 
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For orbital dependent exchange-correlation functionals 
some effort must be put into ensuring that the exchange- 
correlation energy is evaluated correctly. This should 
however be quite straightforward, since all the occupied 
orbitals are known. 



C. Gradients 

Gradients of PESs are easily evaluated in ordinary 
DFT due to the Hellman-Feynman theorem. The 
Hellman-Feynman theorem, however, only applies to 
eigenstates and not linear expansions of eigenstates. Due 
to this there is no easy way of gaining the gradients in a 
linear expansion ASCF calculation. In section IIVCI we 
will show that the Hellman-Feynman gradients do in fact 
not match the true gradients. 



D. Other excitations 

Above we only considered excitations where an elec- 
tron is removed from the Fermi energy and placed in 
some specified orbital. The method is, however, easily 
extended to other types of excitations by representing 
each removed and each added electron as a linear expan- 
sion of KS orbitals. Eq. ([6]) then gains an extra sum for 
each extra linear expansion. In cases of removed elec- 
trons the sign should of course be negative and the sum 
be over KS states below the Fermi energy. Similarly Eq. 
([^ gains extra sums. 



The energy calculation, which is performed after the 
KS equations have reached self consistency, is not sig- 
nificantly different in the linear expansion ASCF scheme 
compared to ordinary DFT. The Hartree energy is eval- 
uated directly from the density, which is also the case for 
the exchange-correlation energy if an orbital independent 
functional is used. So in linear expansion ASCF these 
terms are evaluated exactly as in ordinary DFT. In ordi- 



E. Implementation 

We have implemented the method in gpawf22i^ which 
is a real-space DFT code that uses the projector- 
augmented waves^Si^ (PAW) formalism to represent the 
core electrons. The self-consistent electron density is de- 
termined by an iterative diagonalization of the KS Hamil- 
tonian and Pulay mixing of the resulting densityi^ For 



calculations on single molecules we use the local den- 
sity approximation (LDA)3- as well as RPBE'^^ to de- 
scribe exchange and correlation effects. The LDA is used 
because we compare to TDDFT results obtained using 
the ALDA approximation. ■^'^ and RPBE is used to see 
whether or not the generalized gradient description im- 
proves results. For calculations on molecules at surfaces 
we only use RPBE, because this is designed to perform 
well for molecules adsorbed on transition metal surfaces. 

The projection step described in section Hi Al can easily 
be approximated within the PAW formalism if the atomic 
orbitals are chosen as partial waves; see appendix n for 
details. 

For reasons of comparison we have also made a few 
linear response TDDFT (hTDDFT) calculations. These 
have been made using the Octopus code^i^ which is 
a real-space TDDFT code using norm-conserving pseu- 
dopotentials to represent core electrons. 



III. MOLECULES ON SURFACES 

The linear expansion ASCF method is especially rel- 
evant for molecules on metallic surfaces, because the 
molecular state, due to hybridization, is spread over 
many KS states, ie. it is necessary to write the resonant 
state as a linear combination of KS states. In this section 
we will make a detailed investigation of the 27r resonance 
of N2 on a ruthenium (0001) surface. Furthermore we ap- 
ply the proposed method to several diatomic molecules 
on different metallic surfaces and compare the results to 
other methods and experiments. Finally we map out a 
part of the PESs for N2 on ruthenium (0001) and use it 
to estimate how much energy could possible be put into 
molecular motion from an electron hitting the resonance. 



A. 27r resonance energy for N2 on ruthenium 

The two top panels on Fig. [1] shows the 27r resonance 
energy for N2 on a ruthenium (0001) surface as a function 
of the system size, ie. the surface unit cell and the num- 
ber of ruthenium layers. The resonance energy is the to- 
tal energy difference between a resonant calculation and 
a ground state calculation, both performed with atomic 
positions corresponding to the minimum of the ground 
state PES, ie. it is vertical resonance energies. We mini- 
mize the energy in the ground state calculations by keep- 
ing all surface atoms frozen and found that the nitrogen 
molecule is placed on-top with the two nitrogen atoms 
placed 2.084 A and 3.201 A above the surface. In the 
resonance calculation the 27ry orbital of the N2 molecule 
has been expanded on all KS states above the Fermi en- 
ergy. This expansion has been used as ip^'^^ in Eq. ([S]). 
Although an extra electron is placed on the molecule we 
keep the total number of electrons unchanged, such that 
the unit cell is neutral. This is reasonable because a 




(1,1) (2,1) (2,2) (3,2) (4,2) 2 3 4 5 
Surface unit cell No. of layers 

FIG. 1: (color online) Upper row: The 2tt resonance energy 
of N2 molecule on a ruthenium surface. Lower row: The 
extra charge on the N2 molecule in the resonance compared 
to a ground state calculation. Left panels are for two layers 
and different surface cells, ie. different N2 coverages. Right 
panels are for a (2,1) surface cell and different number of 
layers. The extra amount of charge is estimated using Bader 
decomposition.— i^ 



charged molecule will form an image charge in the sur- 
face, keeping the entire system neutral. 

The resonance energy is converged to within 0.1 eV at 
a surface unit cell of (2,2). The rather large variation in 
energy for smaller unit cells is probably due to dipole in- 
teractions between periodic images. This is confirmed by 
a simple estimation of the dipol-dipol interaction ener- 
gies. The resonance energy is not influenced significantly 
by the number of layers in the ruthenium, indicating that 
the charge redistribution only occurs very near to the sur- 
face. That the charge redistribution is local is confirmed 
by Fig. [21 which shows the change in charge between the 
resonance calculation and the ground state calculation 
for 4 different surface unit cells. For the larger unit cells, 
where the resonance energy has converged, a clearly lo- 
calized image charge is seen below the nitrogen molecule 
and above the first layer of ruthenium atoms. The area 
with extra charge clearly resembles the 27r orbital of ni- 
trogen, indicating that the 27r orbital is well represented 
by the linear expansion of KS orbitals. Fig.[5]also reveals 
that some charge is redistributed within the molecule. 

In order to get an estimate of the size of the charge re- 
distribution we also performed Bader decomposition^i^^ 
on the density found in the ground state calculation and 
the resonance calculation. The two bottom panels on 
Fig. [1] show the extra charge assigned to the nitrogen 
molecule in the resonance calculation compared to the 
ground state calculation as a function of system size. The 
converged value is close to 0.5 electron charge, ie. only 
half of the electron is placed on the nitrogen molecule 
according to the Bader decomposition. This discrepancy 
could either be due to the ambiguity in the way one 




FIG. 2: (color) The change in charge distribution due to the 
excitation. Green: more charge (0.01 a.u. contour), red: less 
charge (-0.01 a.u. contour). The 4 figures are for 4 different 
surface unit cells: (1,1), (2,1), (2,2) and (4,2). Grey atoms are 
ruthenium and blue atoms are nitrogen. The periodic images 
of the atoms are also shown, whereas the density changes are 
only shown in one unit cell. 
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chooses to assign charge to the atoms or a more physical 
effect of charge going back into the surface when extra 
charge is placed on the molecule. The former reason is 
very Hkely, since the image charge is located very close 
to the molecule. 

In order to investigate the effect of charge going back 
into the surface we start by considering the 27r orbital 
itself. Fig. [3] shows the density of KS states and the pro- 
jected density of states (PDOS) for the 27r orbital for the 
ground state calculation and the resonance calculation. 
In the ground state calculation a part of the long tail 
of the PDOS goes below the Fermi energy, ie. a small 
part of the 2-k orbital is occupied here. In the resonance 
calculation the PDOS has moved upward in energy such 
that the tail no longer goes below the Fermi energy, ie. 
some charge goes back into the surface as charge is placed 
on the molecule. Similar effects are seen for the other 
molecular orbitals as visualized on Fig. [4l which shows 
the PDOS for the 3cr, 4(7, Itt and Scr orbitals. Again it 
is seen that all the PDOSs are shifted up in energy as 
more charge is placed on the molecule. Almost the entire 
PDOSs are still under the Fermi level, but small ripples 
can be seen above the Fermi level, also contributing to 
the amount of charge going back into the surface. 

This backtransfer of charge is not an unwanted effect, 
since we try to model the long-lived resonance state, ie. 
the resonably localized peak in the PDOS on Fig. [31 The 
backtransfer of charge is due to some on the energy scale 
very dclocalized bands, indicating a much shorter life- 
time, ie. the backtransfer is expected to happen on a 



FIG. 3: The density of states for a N2 molecule on a ruthe- 
nium slab and the projected density of states on the 27r orbital 
of the N2 molecule. Top: Ground state calculation. Bottom: 
Resonance calculation. 



much shorter time scale than the decay of the resonance. 
It is however clear from Figs. [3] and S] that the charge 
backtransfer in this case is far from the 0.5 electron in- 
dicated by the Bader decomposition. We then conclude 
that the main part of the discrepancy in this situation 
can be assigned to the ambiguity in the way charge is 
assigned to the different atoms. We also find that one 
gets significant different results by assigning charge in a 
different manner, for example by dividing the charge by a 
flat plane midway between the surface and the molecule. 



B. Comparison with inverse-photoemission 
spectroscopy experiments 

In Table |T] we have tested the linear expan- 
sion ASCF method against inverse-photoemission spec- 
troscopy (IPES) measurements and compared the results 
to spatially constrained DFT and ordinary ASCF calcu- 
lations. The modified ASCF values are all calculated in 
exactly the same manner as for N2 on ruthenium in the 
previous section. In all cases the molecules sit on-top, 
and all surface atoms were kept fixed during the mini- 
mization of the molecular degrees of freedom. For the Ni 
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FIG. 4: (color online) Projected density of states (PDOS) on 
the 3(7, 4(7, Itf and 5a orbitals of a N2 molecule sitting on a 
ruthenium slab. The PDOSs are plotted for both the ground 
state calculation and the resonance calculation. The grey area 
indicates energies below the Fermi level. 



TABLE I; Comparison of the 2ir resonance energies for dif- 
ferent diatomic molecules on different surfaces found by spa- 
tially constrained DFT, ordinary ASCF, our modified ASCF 
and experiments. The experimental results have been ob- 
tained from inverse-photoemission spectroscopy (IPES) mea- 
surements: "Johnson and Hulbert,— ''Reimer et al,— '^Reimer 
et al,— ''Rogozik and Dose^. All energies are in eV. We have 
not included IrTDDFT calculations, since it is not applicable 
to periodic systems. 



System 


Constrained ASCF 


ASCF 


Experiment 




DFT 


(orig.) (this work) 




N2 on Ni(OOl) 


2.2 


3.5 


4.0 


4.4" 


CO on Ni(OOl) 


2.2 


3.2 


4.2 


4.Q''/4.5'' 


NO on Ni(OOl) 


2.2 


0.6 


1.4 


i.e^/i.s" 


CO on Ni(lll) 


2.8 


4.3 


4.4 


4.4" 


NO on Ni(lll) 


2.7 


0.5 


1.4 


1.5" 


CO on Pd(lll) 


4.6 


4.1 


4.9 


4.7^ 


CO on Pd step 


2.8 


3.2 


4.5 


4.0'* 



(001) surface we used 3 atomic layers, for the Ni (111) 
and Pd surfaces we used 2 atomic layers. The positions 
of the molecules in their minimized position is given in 
tablelU All resonance energies are vertical from the min- 
imum of the ground state PES. The relevant resonance 
for all the considered systems is the 27r resonance. 

The spatially constrained DFT method was suggested 
by Wu and van Voorhis.— i^^ In the calculations we per- 
form here we divide the space into two areas divided 
by the flat plane mid between the surface and the low- 
est atom in the molecule. We the apply a potential, 
F = I/q . (1 -H exp(^))-i, with Az = 0.2 A and zq 
being the z-value of the dividing plane. Vo is varied un- 
til an extra electron is placed on the molecules side of 



TABLE IL The positions of the molecules in the systems from 
table U] All positions are with relative to the closest surface 
atom. The z-direction is normal to the surface. At the Pd 
step the CO molecule is tilted over the step, which is the 
reason for the composant in the y-direction. All number are 
i Angstrom 

Surface Molecule Pos. of 1. atoms Pos. of 2. atom 



Ni(OOl) N2 
CO 
NO 



Ni(lll) CO 
NO 



Pd(lll) CO 
Pdstep CO 



(0,0,1.638) 
(0,0,1.456) 
(0,0,1.404) 



(0,0,1.774) 
(0,0,1.758) 



(0,0,1.904) O 

(0,0.586,1.801) O 



(0,0,2.798) 
(0,0,2.621) 
(0,0,2.580) 



(0,0,2.941) 
(0,0,2.935) 



(0,0,3.064) 
(0,0.844,2.934) 



the dividing plane compared to the unconstrained cal- 
culation. The energy is then calculated as described by 
Wu and van VoorhiS i^^'^^ The results using the original 
ASCF method have all been obtained by forcing an elec- 
tron in the KS orbital with the largest overlap with the 
27r orbital. 

The results obtained with our proposed modification 
of the ASCF method is seen to agree quite well with 
the experimental results, better than the spatially con- 
strained DFT and the original ASCF method. All the 
results obtained by the original ASCF approach lie too 
low, which is due to the fact that the large hybridization 
of the molecular orbitals makes it impossible to place 
sufficient charge on the molecule. However, a significant 
problem with this method is that PESs often become dis- 
continuous if one chooses to occupy the KS orbital with 
the largest overlap with the molecular orbital, since this 
can be different orbitals at different configurations. 

The major problem with the spatially constrained 
DFT method seems to be that it in some cases is a too 
strict criteria to force an extra electron on the molecule, 
which reflects itself in similar resonance energies for CO 
and NO. We find that the backtransfer of charge dis- 
cussed in the last section is significant for adsorbed NO 
and essential to obtain the resonance energies we find 
with the modified ASCF method. This indicates that 
the spatially constrained DFT approach is more suited 
for systems with a smaller coupling than one has on the 
metallic surfaces considered here. The good agreement 
between our modified ASCF method and experiments in- 
dicate that this method is preferable for these kinds of 
systems and that the backtransfer effect is indeed physi- 
cally reasonable. 



C. Potential energy surfaces for N2 on ruthenium 

On Fig. [5] we have mapped out a part of the poten- 
tial energy surfaces (PES) for a nitrogen molecule on a 
ruthenium (0001) surface in the ground state and the 2TTy 
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FIG. 5: (color) Potential energy surfaces (PES) for a nitrogen 
molecule on a close-packed ruthenium surface in the ground 
state and the 2iiy resonance as a function of the distance 
between the two nitrogen atoms and the distance from the 
surface to the center of mass of the nitrogen molecule. The 
energies are in eV. The small dots represents the points where 
the energy has been calculated in order to generate the sur- 
faces. The black arrow represent a possible trajectory of the 
system in the resonance state (see text). 



resonance. We limit ourselves to two dimensions, which 
at least is reasonable in the ground state, since here it 
is well known that the molecule sits vertically on an on- 
top site. In the resonance state we have tried to rotate 
the molecule a small angle around the surface atom in 
the X- and y-direction at several points on the PES. In 
all cases this leads to an increase in energy, ie. it also 
seems reasonable to stay within the two dimensions in 
the resonance state. Here we will only apply the PES 
to a simple estimate of the possible energy transfer into 
molecular motion from an electron hitting the resonance. 
For a more detailed analysis it is necessary to include 
other dimensions. 

The ground state PES looks as expected, with a small 
barrier for desorption and a local minimum correspond- 
ing to the adsorption configuration. The resonance PES 
has a shifted minimum, which indicates that an elec- 
tron hitting this resonance could induce molecular mo- 
tion, since a sudden shift between the PESs would leave 
the system far away from the minimum, such that the 
atoms would start to move. The maximum possible en- 
ergy gain assuming classical ion dynamics from a single 
electron hitting the resonance can be roughly estimated 
by following the black arrow on Fig. [5] The system is 
most likely situated at the local minimum of the ground 
state PES when the electron hits the resonance. The 
black arrow shows a possible trajectory of the system in 
the resonance state until the resonance decays and the 
system returns to the ground state PES. The potential 
energy after the electron event in this optimal situation is 
approximately 1.5 eV higher than before the event. This 
is seen to be more than enough to desorb the molecule. A 



more detailed analysis involving calculations of the possi- 
ble vibrational excitations and the probabilities of excit- 
ing them will be the topic of a future publication. Such 
an analysis will have to take all six degrees of freedom of 
the molecule into account. 

The PESs show that the center of mass is shifted away 
from the surface when the resonance is occupied. This 
may seem counter intuitive since the charged molecule 
is attracted to the generated image charge in the sur- 
face. However, the resonance weakens the bond between 
the nitrogen atoms, such that the distance between them 
increases, which shifts the center of mass outwards as 
the lower atom is not free to move closer to the surface. 
This effect is more significant than the decrease in the 
ruthenium-nitrogen distance due to the mentioned im- 
age charge effect. 



IV. SMALL MOLECULES 

In the following we present some small tests performed 
on N2 and CO. These small systems have the advantage 
that they make it possible to compare to more accurate 
linear response time-dependent density functional theory 
(IrTDDFT) calculations. When possible we also com- 
pare to experiments. The only advantage of our modified 
ASCF compared to ordinary ASCF for these molecules 
is the possibility of handling degenerate states without 
getting convergence problems, ie. the following should be 
viewed as a test of the ASCF approach rather than a test 
of our modification. We are especially interested in con- 
firming the ability to predict the shift of the minimum 
when going from the ground state PES to the excited 
state PES, which we in section UlI CI argued is very im- 
portant when considering molecular motion induced by 
an electron hitting a molecular resonance. 



A. Excitation energies 

We have used the linear expansion ASCF in combi- 
nation with the multiplet sum method^ to calculate 
excitation energies for different excitations in the N2 
and CO molecule. The results are presented in table 
mil and IIVI respectively. The 4ct and 5a states are 
both represented by a single KS orbital. The Itt and 
2tt states are both double degenerate, so they are both 
represented as a linear combination of two KS orbitals: 
k) = -^Itt K s, a) +'i^\'^KS,b), where iTr^s.a) and \TrKS,b) 
are the two degenerate KS orbitals. The imaginary unit, 
i, has been included in order to get the correct angular 
momentum of the excited states (H and A). This would 
not be possible using traditional ASCF, where one only 
has the freedom to change occupation numbers of the 
KS states. Due to the rotational symmetry of the den- 
sity found from these states the calculations do not suffer 
from any convergence difficulties. That is not the case if 
one just occupies one of the degenerate KS orbitals. Only 



TABLE III: Vertical excitation energies for the N2 molecule 
taken from the minimum energy configuration of the ground 
state. "KS eigenvalue differences, ''linear response calcu- 
lations taken from Grabo et al.'^'*, '^Computed by Odder- 
shede et al^ using the spectroscopic constants of Huber and 
Herzberg.— All theoretical results are obtained using LDA as 
the xc-potential (and ALDA for the xc-kernel in the IrTDDFT 
calculations). 



State Transition Aexs ° 


TDDFT'' ASCF ASCF 
(ALDA) (LDA) (RPBE) 


Exp.'^ 


"'^ 5(7 -> 27r 8.16 
singlet-triplet splitting: 


9.23 
7.62 
1.61 


8.75 
7.55 
1.20 


8.58 
7.52 
1.06 


9.31 
8.04 
1.27 


'^ f^ Itt ^ 27r 9.63 
singlet-triplet splitting: 


10.27 
8.91 
1.36 


10.50 
8.94 
1.56 


10.52 
8.79 
1.73 


10.27 
8.88 
1.39 


"^ 4cr ^ 2tt 11.21 
singlet-triplet splitting: 


13.87 

10.44 

3.43 


11.97 

10.37 

1.60 


12.40 

10.61 

1.79 


13.63 

11.19 

2.44 
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FIG. 6: (color online) The energy as a function of bond length 
for the N2 molecule in the ground state and two excited states. 
The black lines correspond to ASCF calculations, the grey 
(online: light blue) lines correspond to linear response cal- 
culations. The linear response calculations have been made 
using Octopus.— 1^ The vertical lines indicate the positions 
of the minima. 



TABLE IV: Vertical excitation energies for the CO molecule 
taken from the minimum energy configuration of the ground 
state. "KS eigenvalue differences, ''linear response calcula- 
tions taken from Gross et al.— , '^Computed by Nielsen et 
al.— All theoretical results are obtained using LDA as the 
xc-potential (and ALDA for the xc-kernel in the IrTDDFT 
calculations) . 



State Transition Ae/fS " 



TDDFT" ASCF ASCF 
(ALDA) (LDA) (RPBE) 



Exp.' 



5(7 ^27r 



6.87 






singlet-triplet splitting: 



8.44 7.84 7.81 8.51 

6.02 6.09 6.02 6.32 

2.42 1.75 1.79 2.19 



Itt -^ 27r 



9.87 



D^A 

Sa 



singlet-triplet splitting: 



10.36 10.82 10.73 10.23 
9.24 9.72 9.55 9.36 

1.12 1.10 1.18 0.87 



4(7 -> 27r 



11.94 






singlet-triplet splitting: 



- 13.15 13.09 
11.43 12.26 12.09 11.55 

- 0.89 1.00 



the A states are included in the Itt -^ 27r transitions in 
Tables Hill and [TVl since the S states cannot be estimated 
by the niultiplet sum method4^ This is, however, not a 
problem for the kinds of systems for which this method 
is intended, such as molecules on surfaces where high- 
accuracy alternatives are still lacking. 

In general the excitation energies found by the linear 
expansion ASCF method look quite good for the low- 
lying excitations. The accuracy is only slightly worse 
than that of IrTDDFT and significantly better than just 
taking KS eigenvalue differences. The singlet triplet split- 
tings are also rather close to the experimental values. 



The method however seems to struggle a bit more in the 
higher lying 4(t — > 27r transitions. This could indicate 
that the method should only be applied to low lying ex- 
citations. Changing the exchange-correlation functional 
from LDA to RPBE docs not affect the accuracy signifi- 
cantly although a small tendency towards better perfor- 
mance is seen for the higher lying excitations. We note, 
however, that the intended application of ASCF do not 
include simple diatomic molecules, where more accurate 
quantum chemical methods are available. 



B. Excited potential energy surfaces 

The shapes of the potential energy surfaces can in some 
cases be more important than the exact height of them, 
ie. a constant error is not so critical. This is for ex- 
ample the case when considering chemistry induced by 
hot electrons.i^^ In order to get an idea of the accuracy 
with which the linear expansion ASCF method repro- 
duces correct shapes of potential energy surfaces we have 
calculated the potential energy surfaces for the ground 
state and two excited states in the N2 molecule. These 
are plotted in Fig.[6]togcther with results from IrTDDFT 
calculations. 

The small difference between the two ground state 
curves are due to the fact that they have been calculated 
with two different codes. Both codes are realspace codes, 
but gpaw uses the PAW formalism to represent the core 
electrons whereas Octopus uses norm-conserving pseu- 
dopotentials. The calculations have been made with the 
same grid spacing and unit cell size, and with the same 
exchange-correlation potential (LDA/ ALDA). 

The shapes of the potential energy surfaces found from 
the two different methods are seen to be very similar. 
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FIG. 7: (color online) The energy as a function of bond length 
for the N2 molecule in the ground state and two excited states. 
The short thick lines indicate the size of the gradients. 



Especially the predicted positions of the minima are seen 
to agree very well. The shifting of the minima towards 
larger bond lengths is also the expected behavior, since 
an electron is moved from an bonding orbital to an anti- 
bonding orbital. When going to bond lengths beyond 2 
A we start having problems with convergence problems 
in the ASCF calculations, since the 27r orbital ceases to 
exist. This is not a problem we have encountered in the 
systems with a molecule on a surface. 

The good agreement between ASCF and IrTDDFT 
probably reflects that ASCF and ignoring the his- 
tory dependence of the exchange-correlation potential in 
TDDFT are related approximations. For example, the 
density obtained in ASCF would be stationary if evolved 
in time with TDDFT. 



molecules near surfaces, where the molecular orbitals hy- 
bridize, such that none of the KS orbitals can be used to 
represent an extra electron placed on the molecule. The 
method has been implemented in gpawSZiS^ and applied 
to several systems. 

From calculated potential energy surfaces of N2 on a 
ruthenium surface we concluded that an electron hitting 
the 27r resonance in this system can induce molecular 
dynamics due to the different positions of the minima of 
the ground state PES and the resonance PES. Through a 
simple analysis we found that one electron can optimally 
place 1.5 eV in the atomic motion, more than enough to 
desorb the molecule. 

We find good agreement between the model and in- 
verse photo-emission experiments for several diatomic 
molecules on different metallic surfaces. For the consid- 
ered systems we find significantly better agreement with 
experiments using the modified ASCF method compared 
to spatially constrained DFT and traditional ASCF. 

Finally we applied the method to N2 and CO in their 
gas phases we found that excitation energies are esti- 
mated with quite good accuracy for the lower lying ex- 
citations, comparable to that of TDDFT. Especially the 
shape of the potential energy surfaces and the positions 
of the minima agree well with TDDFT results. 
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APPENDIX: OVERLAPS USING PAW PSEUDO 
WAVEFUNCTION PROJECTIONS 



As mentioned in section III CI the Hellman-Feynman 
theorem does not apply in the linear expansion ASCF 
method. This is verified by the calculations shown on 
Fig. \7\ Here the energy of the ground state and two ex- 
cited states in the N2 molecule is plotted as a function of 
the bond length. The short thick lines indicate the gra- 
dient given by calculated Hellman-Feynman forces. For 
the ground state the agreement is as expected perfect, 
but for the excited states there is a clear mismatch. Un- 
fortunately this implies that it is computationally heavy 
to do dynamics or minimizations in the excited states. 



The Projector Augmented Wave (PAW)^ method uti- 
lizes that one can transform single-particle wavefunctions 
iV'n) oscillating wildly near the atom core (all-electron 
wavefunctions), into smooth well-behaved wavefunctions 
\ipn) (pseudo wavefunctions) which are identical to the 
all-electron wavefunctions outside some augmentation 
sphere. The idea is to expand the pseudo wavefunction 
inside the augmentation sphere on a basis of smooth con- 
tinuations {(pf) of partial waves {(pf) centered on atom a. 
The transformation is 



i^„) = i^„) + E(i'^")-i'^°))(p°i^") 



(A.l) 



V. SUMMARY 

We have extended the ASCF method of calculating 
excitation energies by allowing excited electrons to oc- 
cupy linear combinations of KS states instead of just sin- 
gle KS states. This solves the problems encountered for 



where the projector functions |j5") inside the augmenta- 
tion sphere a fulfills 



Y^mm^h {p-\4>-) = s.„ 



|r-R°| <r^. 
(A.2) 
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Suppose we have an atom adsorbed on a metal sur- 
face and we wish to perform a ASCF calculation where 
a certain atomic orbital \a) is kept occupied during the 
calculation. If the orbital is hybridized with the metal 
states we need to find the linear combination which con- 
stitutes the orbital. This can always be done if a suf- 
ficient number of unoccupied KS orbitals is included in 
the calculation 



\i) = ^C„i|V'n), Cni = {lpn\ 



(A.3) 



Since the partial waves are typically chosen as atomic 
orbitals we just need to consider the quantity 

{^^m = {^nm +J2{^n\pf){{ci^fm-{4>f\<i>r> 



j,a' 



If we were just considering a single atom, the last equality 
would be exact inside the augmentation sphere since the 
partial waves would then be orthogonal and the pseudo 
partial waves are dual to the projectors in Eq. (|A.2p . 
When more than one atom is present there is correc- 
tions due to overlap of partial waves from neighboring 
atoms and non-completeness of projectors/pseudo par- 
tial waves between the augmentation spheres. However 
using {ipn |p°) is a quick and efficient way of obtaining the 
linear combination, since these quantities are calculated 
in each step of the self-consistence cycle anyway. The 
method can then be extended to molecular orbitals by 
taking the relevant linear combinations of ('0„|p°). 



ii'nlpt)- 



(A.4) 
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